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Abstract 

The three-loop QCD corrections to the p parameter from top and bottom quark 
loops are calculated. The result differs from the one recently calculated by Avdeev 
et al. As function of the pole mass the numerical value is given by Sp = (1 — 

2.8599 ^ - 14.594 (^) 2 ). 



1 Introduction 



The precision of electroweak observables measured at LEP, SLC and the TEVATRON has 
stimulated a variety of theoretical calculations. These are required to match the exper- 
imental accuracy and to pin down the parameters of the Standard Model, in particular 
the top and the Higgs mass and to search for new physics. 

A cornerstone in this analysis is the evaluation of top quark contributions to the p 
parameter. With the high value of M t as suggested by the CDF-collaboration |1| and the 
strong sensitivity of p to small variations of M t through the quadratic dependence 
precise theoretical predictions become mandatory. 

Top mass corrections to the Zbb vertex, the only other place with a strong dependence 
on M t ||, |4]] , are specific to this reaction, while those to the p parameter enter numerous 
relations between observables. This is the second justification for a precision calculation. 

In addition to the 0(GFM^a s ) two loop contribution || (for a related calculation 
based on dispersion relations see ||, for a recent review see 0), the two- loop electroweak 
corrections have also been evaluated, for vanishing M and even for arbitrary Mjj 0- 

In this paper the three-loop result of order GpM^a 2 s is presented. A similar calculation 



has been performed by Avdeev, Fleischer, Mikhailov and Tarasov ||10|| . However, our 
calculation disagrees with their formula. In the following section 2 details of the calculation 
will be described, in section 3 the result will be presented and a brief numerical discussion 
given. 



2 The calculation 

Quantum corrections to the p parameter can be connected to the gauge boson self-energies 
through 

n z (o) nx'(o) 
''"-ST" ST" () 

Here n(0) denotes the transverse part of the polarisation tensor I\.^ v {q) for vanishing 
momentum q. The evaluation of these self-energy diagrams is performed for = and 
within the framework of dimensional regularisation. Large intermediate expressions are 



treated with the help of FORM 2.0 For D ^4 anticommuting 75 was used, except 



for the double triangle diagram. In order to evaluate this diagram, which is related to the 



axial anomaly, the definition of 't Hooft and Veltman |fL2 |, formalized in [[y|, was applied. 
Its contribution is finite and the result coincides with a previous calculation |14| where 
D = 4 from the very beginning. A covariant gauge with arbitrary gauge parameter for 
the gluon propagator was chosen. 

The tadpole integrals required to calculate the one- or two-loop corrections are easy 
to evaluate even for arbitrary powers of the propagators. This does not hold true for the 
three-loop case. After performing the traces the reduction to scalar integrals is performed 
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by decomposing the scalar products of the numerator in appropriate combinations of the 
denominator. Subsequently, recurrence relations provided by the integration-by-parts (IP) 
method [[DJ are used in order to reduce every scalar Feynman integral to a small number 
of so-called master diagrams which have to be calculated explicitly. The IP method was 
first applied in [IB| to three-loop tadpole integrals. There the subclass of those diagrams 
which contain a continuous massive quark line and which are relevant for the Z boson 
self-energy was considered. For the calculation of the p parameter the method has to be 
extended to a second class of integrals originating from the W self-energies. One thus 
arrives at three master integrals 
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(The same variables and conventions as those of [0 are adopted. Following standard MS 
practice we discarded terms proportional to 7^ and log47r in the r.h.s.) The first master 
integral has been calculated analytically in [0; the results for the last two integrals 
are given |ITJ (D 3 is presently only known numerically). We checked the result for 
numerically. For 5*2 we reproduced the analytical result. The evaluation of D3 is described 
in [|T7| . Employing a different method, we obtain the result given below, which is consistent 
with |fr7| . The values for the constants -B4, S2 and D 3 follows: 
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3 Results and Discussion 



After performing mass and charge renormalization in the MS scheme the following result 
for the W boson propagator is obtained: 
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The coefficient of Cp((4) differs from the recent result of [[K]]. Its value 26 has to be 
compared with 88/5. This leads to a significant modification of the numerical predictions 
to be discussed below. In this expression rif = 6 denotes the total number of quark 
species and Z = log /i 2 /m 2 . The result is expressed in terms of the MS renormalized top 
mass m t (/z 2 ). The variable x t is defined as 
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From now on the explicit //-dependence both in x t and a s will be suppressed. From the 
context it should be evident which scale is adopted. 

For the Z boson propagator the following result is obtained 
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in agreement with [fLOR . Both II^ and II z are independent of the gauge parameter. As 
expected this holds true even before mass and charge renormalization are performed. 

Eqs. (@) and @ immediately lead to Sp in terms of the MS renormalized top mass. 
The poles which are still present in the Z and W self-energies individually cancel. 
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With the help of the relation between the OS and running top mass [18|] for p? 
m t {M 2 t ) = M-i 
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the result is easily expressed in terms of the OS mass: 
$Pos = 3X t 
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Here M t is the pole mass and X t = GfM 2 /8^/2ir 2 . The residual log p terms are cancelled 
by the /z-dependence of a s . 

At this point two consistency checks should be mentioned which were performed in 
order to test the correctness of our result. The first one is a different method of calculating 
bp respectively the polarisation functions for the W and Z boson. It relies on the axial 
Ward identity which connects the axial part of the polarisation tensor with the pseu- 
doscalar polarisation function. (The double triangle diagram was not considered in this 
context.) If the fermions in the loop have the masses mi and mi the following identity 
holds: 



The second term of the r.h.s. is independent of q and therefore not relevant in this context. 
The l.h.s. in lowest order is already 0(q 2 ). IP(g) was evaluated up to 0(q 2 ). Because of 
the different tensor structure, different recurrence relations have to be applied to compute 
the p parameter. IP has to be expanded in q up to order q 2 . Additional propagators are 
generated and different parts of the FORM programs become relevant. In addition it was 
checked that the q independent part of IP cancels against the second term such that the 
r.h.s. indeed is of 0(q 2 ). 5p as given in eq. fl5|) was reproduced with this method. 

The second check is also connected with an expansion in the external momentum. 
The polarisation tensor of the vector bosons was expanded up to order q 2 . The external 
momentum was routed through the graphs in two different ways. Again the same result 
was obtained for every diagram although the intermediate steps are very different. 

Substituting Ca = 3, Cp = 4/3 and p? = m 2 or p 2 = M 2 in eqs. (||) and (0) respec- 
tively a fairly compact form for 5p is obtained: 



q^q v W v ' a (q) = (m 1 + m 2 ) 2 W(q) + (mi + m 2 ) < 0|gigi + q 2 q 2 \0 > . 
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To evaluate these results numerically the values for B4, S2, -D3 and the ( functions are 
inserted. For convenience of the reader we display separately four different contributions: 

(i) the contribution from the double triangle diagrams related to the axial anomaly JTj |, 

(ii) the contributions with exactly one fermion loop which together with the previous one 
would give the result in the "quenched" approximation, (iii) the contribution from light 
quark loops proportional to rii = n/ — 1, and finally (iv), the contribution with the light 
quarks replaced by top quarks. 
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The coefficient in front of ni in eq. ( |T2D is in reasonable agreement with the numerical 
result in |19|]. Eqs. ( p]) and ( |i~2"D in particular the separation of the various contributions 
allow the use or test of a variety of optimization schemes which is left to the reader. 
For the final result and after setting rii = 5 we obtain 



6 Pms = 3x t \1 - 0.19325^ - 3.9696 
5p os = 3X t [ 1 - 2.8599 — - 14.594 (— 
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(14) 



The coefficient in front of a 2 s in the MS result differs significantly from [10|: —3.969 . . . ver- 
sus +0.07111 . . .. The OS results differ by the same amount —14.59 . . . versus —10.55 . . .. 



In JTIJ it was mentioned that the contribution from the double triangle diagram, 
associated with the axial anomaly alone amounts to about 40% of the total three- 
loop correction. Here a fraction of approximately 30% is still traceable back to this single 
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diagram in the OS scheme. In the MS scheme the corrections are completely dominated 
by this diagram. 

Finally the numerical effect on the prediction for My/ and sin 2 e g- from a,G F and 
Mz will be discussed. If sub- leading terms are neglected the following relations hold [2T| : 



M 2 W 



sin 2 cff 



pM 2 




Here A = ^Txaj\f2G F = 37.2802 GeV, Aa » 0.06 and p = 1/(1 - dp). The relative size 
of the one-, two- and three- loop corrections is given in Table [l]. The first row indicates 
the size of 5p itself. Rows two and three give the relative contribution of the one-, two- 
and three-loop corrections with respect to the Born result in the OS and MS scheme. The 
numbers are obtained with the following input data: a s (M 2 ) = 0.1092, m t (M 2 ) = 164 
GeV, M t = 174 GeV, M z = 91.188 GeV and G F = 1.16639 10" 5 GeV" 2 . The value of 
a s (M 2 ) = o4 6) (M 2 ) was obtained from af\hl 2 z ) = 0.120 with the help of the three-loop 
(5 function and the matching condition between af\p?) and af\p?) at the scale 
p? = M 2 [2^|. (The numbers in brackets indicate thereby the numbers of active flavours.) 
The MS top mass mt(M 2 ) was derived from eq. (§). This value serves as starting point 
in order to calculate m t (p 2 ) using the corresponding three-loop evolution equation |23| . 
(Unlike the coupling constant the running top quark mass is only defined in full n/ = 6 
theory.) 



OS 


+l-loop 


+2-loop 


+3-loop 


5p 

5M W /M W 

5 sin 2 O e g/ sin 2 O e g- 


0.00949 
0.00682 
-0.01349 


0.00854 
0.00614 
-0.01216 


0.00838 
0.00601 
-0.01192 




MS 


+l-loop 


+2-loop 


+3-loop 


Sp 

5M W /M W 

5 sin 2 O e ff/ sin 2 O c g- 


0.00843 
0.00605 
-0.01200 


0.00837 
0.00601 
-0.01191 


0.00833 
0.00598 
-0.01186 



Table 1: Numerical results including successively higher orders. 

It is interesting to mention that the scheme dependence decreases enormously when 
taking higher loop corrections successively into account. 
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One observes that the three-loop corrections are (at least for M\y) not at all small 
compared with the two-loop QCD contribution. Furthermore they are approximately of 
the same order of magnitude as the two-loop electroweak result because for mnjm t = 1.5 
the contribution from the two-loop electroweak correction amounts —0.018% [T^, § to be 
compared with the three-loop QCD corrections of —0.013%. It is also possible to translate 
the 0{a 2 s ) corrections directly in a change of the top mass contained in X t . In the on-shell 
scheme this corresponds to a change of approximately —1.5 GeV. 

The dependence of the result on the renormalization scale p? is shown in Figs. |l|a and 
[TJb for 5pos and Sp^ respectively. The same input parameters have been used as before. 
The dotted line in Fig. [TJb gives the one-loop prediction (which is constant for Spos and 
completely off-scale in Fig. Qa), dashed and solid lines represent the two- and three- loop 
results. The prediction is clearly stabilized through inclusion of higher orders. 

Another possibility would be to absorb the a 2 s contribution in the choice of an effective 




Figure 1: Renormalization scale dependence of Spos (a) and Spy^ (b). 
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scale of the a s correction: 

/ a s ((0.302M 4 ) 2 )\ , . 

5p os = 3Xt M - 2.8599 — — — J. (15) 

To summarize: The evaluation of the three-loop QCD correction to the p parameter 
has been repeated with a result different from the one of [IIJ . The numerical difference is 
sizeable. 
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